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Abstract 



We consider two different problems involving tlie localization of a single 
polymer chain: (i) a periodic AB copolymer at a selective fluid-fluid inter- 
face, with the upper (resp. lower) fluid attracting A (resp. B) monomers (ii) 
a homopolymer chain attracted to a hard wall (wetting) . Self avoidance is ne- 
glected in both models, which enables us to study their localization transition 
in a grand canonical approach. We recover the results obtained in previous 
studies via transfer matrix methods. Moreover, we calculate in this way the 
loop length distribution functions in the localized phase. Some finite size 
effects are also determined and tested numerically. 
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I. INTRODUCTION 



The behaviour of polymers at interfaces is currently a subject of great activity™. In 
this paper, we focus our interest on two problems, namely (i) the adsorption of a periodic 
copolymer chain on a selective fluid-fluid interface (ii) the adsorption of a homopolymer chain 
on a hard wall (wetting). These problems have been previously solved for one dimensional 
geometries, provided certain assumptions are fulfilled (self avoidance is neglected, wetting 
is considered in the SOS approximation,...). A quick guide to the recent litterature can be 
found in referenced!^ and references therein. Here, we will tackle these models, with the 
same assumptions as above, through the use of a grand canonical method. This method 
was initiated in the context of the DNA helix-coil transitionUll and was recently used 
in the context of RNA foldinglii. We will recover previously known results about models 
(i) and (ii). Moreover, this approach enables us to obtain the loop length distribution 
function in the localized regime. We also find in this way some analytic expressions of finite 
size corrections at the critical temperature, which are in agreement with numerical transfer 
matrix calculations. The lay out of the paper is as follows. Section || deals with the periodic 
AB copolymer chain at a selective interface. In section PTT] , we consider the wetting of a 
homogeneous wall by a homopolymer chain. 



II. PERIODIC COPOLYMER CHAIN AT A SELECTIVE INTERFACE 

A. Introduction 

We consider a periodic AB copolymer chain of 2N monomers at a selective fluid-fluid 
interface, located at = |- The upper (resp. lower) fluid is a good solvent for A- type (resp. 
B-tjpe) monomers, and will herefrom be called fluid a (resp. fluid f3). This situation is 
modelled through the Hamiltonian 

2N 

T-(-2N = - E sgn[z{i) - -] (1) 
i=i ^ 

where sgn('u) denotes the sign function, and the "charges" qi of the periodic chain are chosen 
as g2n+i = Q'A > and q^n = —qs < 0, together with > g^. In equation (|l|), z{i) denotes 
the position of monomer i of the chain, which we represent as an (unidimensional) random 
walk (RW) with the properties z{i) = .., —2, —1, 0, +1, +2, .., and z{i + 1) — z{i) = ±1. We 
emphasize once more that self avoidance is neglected in this model. 
The partition function then reads 

Z2N = E e-^"^^^ (2) 
(RW) 

where (3 = -f and the sum runs over all random walks (RW) starting with an A monomer at 
z{i = !) = !. The thermodynamical properties of this model can be obtained through the 
use of transfer matrix methodJl, as given in Appendix A. Here, we follow referencesiHii, 
and consider a particular chain configuration as a collection of non interacting loops. This 
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class of problems is solved in a grand canonical approach, by defining a grand canonical 
partition function 

+00 

Z{K) = Y: K'^Z^n (3) 

N=l 

and noting that the radius of convergence K^{f3) of the series Z{K) gives the free energy 
per monomer in the thermodynamic limit 

■'2(3 ^ ^ 

Moreover, one can easily calculate in this way quantities such as the loop length distri- 
bution function, which are not accessible through the transfer matrix approach. 



B. Calculation of the grand canonical partition function 

The partition function Z2N for a chain of (2A^) monomers can then be decomposed 
according to the position z{2N) of the end point of the chain. For simplicity, we set z{2N) = 
2n, and write 

+00 

Z2N= E Z2N{2n) (5) 

n=— 00 

Given the expression of the Hamiltonian (|T]), the chain configurations can be further de- 
composed into loops. A preliminary remark is as follows: if A/'(/, z) denotes the number of 
random walks of / steps going from z = to 2; > in the presence of an absorbing barrier 
at 2 = 0~, the image methodo yields 

Ar(2/,2p) = C^+^-C^+^+^ (6) 
AA(2/-l,2p-l)=C^tV-C^r-"i (7) 

depending on the parity of {l,z). 

Let us temporarily suppose that the chain end is in fluid a {2n > 2). The first loop, 
located by hypothesis in fluid a contains (2/i — 1) monomers, namely li A monomers and 
(Zi — 1) B monomers. The second loop, located in fluid (3, contains (2Z2 — 1) monomers, 
namely (/2 — 1) A monomers and I2 B monomers.... The last loop is denoted by l2p- Finally, 
one must also specify the number 21 of monomers between the interface and the chain end. 
We may then write 

+00 00 00 00 00 2p 

^2^(2n)= E E E - E E'^(2A^-(E(2^^-l) + 20)e^''"'"^^''""^"^'W(2/i-2,0) 

p=0«i=li2 = l l2p = ll=l i=l 
^PqA-P{qA-qB)l2J^^2l2 - 2, 0) e'^'^^' '^'^'^^''^^^^^^ M {2l2p - 2, 0) 

e/3(qA-te)/_yv-(2/ - 1, 2n - 1) (8) 
Plugging equation (H) into equations @ and (|^), we get 
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Z{K)= Z{K,2n) (9) 

n=—oo 

where 

+00 

Z{K,2n) = K''Z2N{2n) 

N=l 

-1 00 

with 

9o = 7^ (11) 



and 



+00 



Biy) = E y^'- W(2/ -2,0) = -[l- ^1 - V j (12) 
We finally get, for (2n > 2) 

Z(i^,2n) = X:^ir^Z2A.(2n) = ^ _ e/3(,,+,,)ij(;^i/2e/3,o)5(iri/2e-/3*0 ( iTi/^e/^-^o " ^^3) 
Similarly, if the chain end is in fiuid (3 {2n < 0), we find 



C. Thermodynamics 

We finally get from (|13D and (|T|) 



+00 



Z{K)= y Z{K,2n) = "^^'^^ " ' " , , ,7/; , \ V/^^!, \" ^ (15) 



1 _ (,0(qA+qB)B{K^/^el^'io)B{K^/^e- 



where 



r (y) = ^^^^ = y (17) 

y 

The radius of convergence K^:{(3) of the grand canonical partition function Z{K) is 



either given by the singularity of the numerator of equation ([I5D , or by the singularity of 
the denominator. The former corresponds to the value 
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KaeUfi) = \e-^^'^' (18) 

The latter comes from the resummation of the loops, and is given by 

sinh(/?gA) sinh(/jgB) 

^*KP) = . ,2 Of ; (19) 

smh p{qA + qB)) 

If the solution K^{(3) is smaller than K^eioc then the chain is localized. This actually 
happens for (5 > (5c (i-e. T < Tc), where /3c is the solution of 

g2/3cte ^ g-2/3cgA _ 2 = (20) 

In particular, in the limit go = (q'a — 9b)/2 — > 0, the critical temperature Tc diverges as 

T ~ — (21) 

" 90-0 go ^ ' 

Note that the scaling of Tc with go differs from equation (10) of referenceB, since the model 
used in this reference allows the monomers to sit at the interface, in contrast to equation 
(0). The physics is nevertheless the same: a symmetric copolymer chain (g^ = qs) is always 
localized at the interface. 

Since the radius of convergence K^[(3) and free energy per monomer are linked by equa- 
tion (H), we get 



1 ^JsmhjPqA) sinh(/?gB) 
^^^^ = P'^ sinh/5(g. + g.) ^''^ 

for T < Tc, which in turn yields 

f{T) ^^^^ Uiocin) - C{Tc - Tf (23) 

where fdeiociTc) = —Tc In 2 — go is the free energy of the delocalized phase (since we have 
chosen go > 0, the chain is then in fluid a). Equation (^) corresponds to a critical exponent 
a = 0. The jump in the specific heat per monomer is given by 

^ .^-^ _ ^^ (26^^-"^ - l)ln(2 - e~'^^^-) - 2/3cg^ V 

^-(^^ ) - [ 2(e2/3.<?>. - 1) J ^^^^ 

Coo(T+) = (25) 



D. Density profile 

We denote the normalized density profile by p{2n) and write 

p(2n) = hm (26) 



^2N 
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Inverting equations (|T3D and (0), we have 

/ dKZ{K,2n) 
ZM2n) - f^.^^^^ ^^^^^^ Q ^-^^ ^27) 

After deformation of the contour, the dominant contribution for large comes from the 
pole at K^{(3) (see equation (|T^), and reads 

Z (2n) ~ dKZ{K,2n) ^^g) 

Aflarge /circle around K^, 2m K^+^ 

Summing over n leads to 

Z ^ i dK Z{K) ^29) 

Tvlarge /circle around AT* 2i7r Ar^+^ 

The density profile p{2n) is then obtained as 

1 / B(K^^'^e'^'^°) \" 
P(2^) = 175 ^ T79 T79 ^T79 " 1 (30) 

for n > 1 and 

for n < 0. Inserting the value of K^{[3), we get 

_ (e2/3te + e-^/^'^A _ 2) (e2/35A + e-^/^^s - 2) / 1 - e'^/^^A \ " 
P^'^^> ~ [e^l3{qA+iB) - i)(i _ e-2/39A)(i _ e-2/3<?s) e^/^-Js - 1 j for n > 1 (32) 

^(2^) - (e2/3(g^+gs) _ _ e-2/3gA)(i _ e-2/5<?B) e2/39A _ 1 j for n < U (ddj 

i.e. the density around the interface decays exponentially from the interface with char- 
acteristic lengths (resp. (,13) in fluid a (resp. fluid P), with 

UT) = -T^^ (34) 
UT) = , / X (35) 



Close to Tc, we get 



and 



^a{T) oc ^ ^ ^ (36) 

UTc) = , . (37) 
(^A + ge) 

The order parameter of the transition can be chosen as the probability to be in fluid 
/? and vanishes linearly at the transition: 

M(,= p(2n) ^ocjT, - T) (38) 



n=—co 
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E. Probability distribution of loop lengths 



Contrary to the direct matrix transfer approach, the grand canonical approach allows 
one to obtain informations on the statistical properties of loops in each fluid. Without loss 
of generality, we may consider the partition function Z2Ar(2n = 0) where the chain of length 
2N is attached at the interface at both ends, and decompose it into 

oo oo 

Z2N{2n = ^)=Y.ll Z2N{2n = 0, 2h - 1, 2/2 - 1) (39) 

where Z2Ar(2n = 0, 2/i — 1,2/2 ~ 1) represents the partition function with the constraints 
that the first loop in fluid a contains exactly (2/i — 1) monomers and that the first loop in 
fluid (3 contains exactly (2/2 — 1) monomers. 
We again define generating functions by 



00 00 00 



Z{K, s„ ^2) = E E E i^^e-^^(2'-i)-^^(2'-^)Z2^(2n = 0, 2/1 - 1, 2/2 - 1) (40) 

Ar=l Zi=l l2=l 

The generalization of equation ([T^) to the case (si,S2) 7^ (0,0) then yields 

e/3{gA+te)5(;^i/2g/3goe-«i)5(i^i/2e-/3goe-«2) 

Z{K, Si, S2) = ^ _ 5(;^-l/2g/3go)5(7^1/2e-/3go) V'^^) 

Again, by inversion and deformation of the contour in the complex plane, we get the 
asymptotic behavior for large 



Z2n{Si, S2) = f 

r 7.iK e, 

(42) 



circle around 1m K^+^ 

dKZ{K,suS2] 



Marge /circle around K^.. liir K^^^ 

The generating functions Pq, and Pg of the loop length distribution functions in the two 
fluids are given by 

ij^i ^-^\Z2m{s^ = 0,S2 = 0)J ^ ' 



00 



P,(.) ^ Y. e-'<-'-"m2k - 1) = m„ |Hii|£i^^lil^ (44) 

h=l V^27v(Sl = 0,S2 = 0)/ 

More explicitely, we have 

^(fCy'e^^oe-^) , 1 - Vl - e--('3)-2- 
Pais) = — 7775 = , (45) 

p. s = ^ _ L = e' , , (46 
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with 



,^2l3{qA+qB) 



.{(3) = -2/?go - 21n2 - InK^ = 2/?(/,,„, - /) = In l^ ,^,,,;^^,,,, _ _ ) (47) 

We note that equations ( ^Sf ) and (|46|) exhibit several loop size scales, that we consider 
below. 

At the critical temperature (T = Tc), the distribution of loop lengths in fluid a is 
simply the loop distribution of a free random walk 



Pais) =e' 

T=Tc 



1 - Vl - = 1 - V2s + 0(s) (48) 
yielding an algebraic decay for large 1^ 



(■a 



A consequence of this critical distribution (T = Tc) is as follows. For a chain of length N, 
the number of loops scales as V^, the longest (a) loop being of order N. 

In the critical region T ^ where a;(/3) ~ (T^ — T)^, the moments of (/„) diverge as 

oo / 1 

/S = E(2^ - l)"^"(2/ - 1) - T) [jjr-^,j (50) 

i.e. for large the correct scaling variable is A = la{Tc — T)'^, but the normalisation of the 
region of flnite A varies as {Tc — T). This can be understood by considering the normalisation 
of the region of loops bigger than a given scale Ic ~ ^^.^2 for a free random walk ( ^91) 

oo 1 1 

(51) 

In other words, the phase transition is driven by a small fraction of large loops. 
It is of interest to note that the weight of the small loops is flnite, even at Tc. For 
instance, the probability to have /q, = 1 is 

pn = 1) = — -— - - (52) 

At zero temperature we of course recover the ground state 

Pa{L)^^^5l^,l (53) 

In marked contrast, the critical distribution of loop lengths in the (jS) fluid is given by 



1 — J] — p~4:l3cqo-2s 

Pa{s) ~ e'- ^ (54) 



leading to flnite moments 
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F. Correlation function 



To further emphasize the meaning of the scale corresponding to uj{f3) in equation 
we consider the connected correlation function 

11 11 

C{j - i) =< sgn[z{{) - -] sgn[z{j) --]>-< sgn[z{i) - -] >< sgn[z{j) - ^] > (55) 

Its generating function reads 

oo 

C{K)=J2K^-'C{j-t) (56) 
j=i 

and can be expressed, via a loop decomposition, in terms of the functions Pa{s) and P/sis) 
defined in (|^,I^), with the replacement [s = —InK). Inverting (|56|), we have 



CU !^4i^ (57) 

Jcircle around 2in K^^^ * 

At large (j — i) separation, the behaviour of C{j — i) is dominated by the singularity at 



K = e 2 of Pq,(— Ini^), leading to 



C(j-z)^ oc e-(^-^)^ (58) 
This shows that uj{i3) is the inverse of the correlation length along the chain. 

G. Finite size properties 

In this section, we again consider that both chain ends are fixed at the interface. This 
means in particular that one has Z2N = ■Z'2Ar(0), Z{K) = Z{K, 0), 

1. Free energy 

In the inversion formula of equation (pT]), we set 2n = and get 

Z (0) - </ dKZ{K, 0) ^^g^ 

/circle around 2i7r K^~^^ 

where Z{K,0) given in ([l^) presents a simple pole at K^{P), and a cut on the real axis, 
namely [KdeiodP), +00]. 

In the localized phase, and after deformation of the contour, the dominant contribu- 
tion for large N comes from the pole at K^{i3), and the correction is given by the leading 
contribution of the cuts. Simplifying as explained above the notations, we have 

'j^m.o(^] (60) 



Marge /circle around i^* 2m K^^^ 



deloc / 



' ^P{0)+0(-^Y] (61) 



Marge i^fV 



deloc 
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where p(0) given in (^) represents the probabihty to be at 2n = 0. 

In the delocahzed phase, the finite size properties are given by the contributions of the 
cuts 



r+co f^jr X 

Jo 2m I (l + xF+i [Z{K,,ioc{l + x + to)) - Z{K,Ml + x- to))] (62) 



To obtain the dominant behavior at large A^, we will set v = A^ln(l + x). 
For T > Tc, using 

Z{K,,,U^ + X + lo)) - Z{K,,U^ + x- lo)) = (1 + 0{x)) (63) 

with 

G{[5) = e^l^^^ (1 - .^1 -e-4/39o) _ 1 (64) 

we get the asymptotic behavior 

'■+°° dv 1 
'^^jvlarge ■^'o 2A^z7r K^^i^^ 

^ 1 (l + g(/3)) 

ivlarge K^iocV^N^^ 2G(/?)2 

For T = Tc, G{(3c) vanishes and we have 

Z{Kdeloc{l + X + lo)) - Z{Kdeloc{l + X - lo)) = {I + 0{x)) (67) 

leading to 

'■+°° dv 1 
' Marge -/o 2Nm K^^^^^ 



/ TTTT-T^e-^ Z(ir,eZoc(e^+^o))-Z(ir,eZoc(e^-zo)) (65) 

(66) 



^2A/ ^ / " ^J. e-^ [z(ir,e«oc(e^ + ^o)) - Z{Kdeioc{e^ - ^o))] (68) 



1 



jvlarge K^^i^^V^ 



(69) 



This shows that the probability of presence at the interface decays with the length as 

at T = Tc (H), whereas it decays as l/A^3/2 ^ j.^ 

In the following, we will need the general scaling form which contains (|66|) and (|69D as 
special cases 

..i+^fL, , f v^ml (70) 



Marge /<-iJ,„yA'ff(/3) ' 

where the scaling function q{x) reads 



g X = - / dv—^ 71 

TT Jo f + X"' 
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and 



H{(3) = ^ , ^ (72) 



and where has been defined in (|6^). 

In particular, near T^, we have for T > the scahng form 

leading to the finite size free-energy 



5. Specific heat 

In the thermodynamic limit, the specific heat Coo(T) = —T '^J^f^ presents a jump at 
given in equation (^). To compute the finite-size specific heat C2Af(T'c) at T = Tc, we use 
the scaling form ( [74D 



^ ^^^^2^1 _iMo) /g'(0)y\^ , g'(0) ^/^(^ 

C^2A.(Te)--T^|^=^,-- -^-(^-^j K^oo(Tj + ^ ^ +.. (75) 



-^C^{T^ ) - Vtt ^ — + .. (76) 



Similar relations between Con{Tc) and 6*00(71;) are found in phase transitions displaying 
a jump in the specific heatoii^EJ. 



linked 



3. Order parameter and uniform susceptibility 

A possible order parameter at the transition has been mentionned in equation 
ied to the probability to be in fiuid (3. We consider here 



2N 
M2N = 

i=l 



ilD, and 



(77) 



and the order parameter m2N = ^^^f^- Note that m2N is related to the quan 
(^) by {m2N = 1 — 27V1/3). Given the expression of the Hamiltonian (|lD, we get 



moo =< sgn 



zii) - ^ >= -dhfiT, qa^ qA + h, Qb ^ Qb - h)\h=o 

^ sinh^(g^-g,) ^ ,_ctein_T) + 
2 smh pqa smh (JqB t^t- 



(78) 
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In the thermodynamic hmit, the corresponding susceptibihty reads 



I- (< > -(< M2N >)') = -dlf{T, qA-^qA + h, Qb ^ Qb - h)\h=o 



= - f ' + ' I (79) 

2T \^ sinh^ pq^ sinh^ Pq^J ^ ^ 

leading, at criticahty, to 

= (1 _ e-l..)2 (80) 

To compute the finite-size behavior of the order parameter and of the susceptibihty 
X2N, "we use again the scahng form ([7^), but where now Tc is a function of h defined by the 
equation for Tc where qA ^ qA + h, qB ^ qs ~ h. Since the calculations exactly parallel the 
ones described above, we only quote the results 



and 



, , V27r 1 , , 



X2n{Tc) = ^Xoo(T-) (82) 



4. Numerical results 

We have done numerical calculations using transfer matrix methods with both chain 
ends fixed at the interface. We chose g^i = 1 and = 0.5, yielding a critical temperature 
Tc = l/(ln(i^i^)) = 2.07809... (see equation (|20|)). The corresponding values of the critical 
specific heat and susceptibility are CcxiiT^) = 0.11056... and Xoo(T~) = 5.03932... (see 
equations (H), (|0D). 

Results for C2n{T) and X2n{T) in the critical region are shown in Figures 1 and 2, for 
chain length up to 2N = 18000 and are in quantitative agreement with the above values. 
We get for instance a value C2n{Tc) ~ 0.04739.. as compared to the value 0.047454... coming 
from equation ([T^), and X2NiTc) ~ 2.13376... as compared to the value 2.16289... coming 
from equation (^2]). We have also tested equation (|8TD : Figure 3 shows the behaviour of 
-'m2N{T)) critical region. At Tc, the theoretical value coming from equation 




2 

(^) is 2.0279... as compared to the "experimental" value 2.0130.. of Figure 3. 
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III. THE WETTING TRANSITION 



A. Calculation of the grand canonical partition function 

We consider the SOS version of the 2D wetting problem near a attracting wall at h = 0. 
Our presentation will be rather brief, and we follow the notations ofi andi. Let us denote 
by {h{i),i = 1,2,..,A^} the consecutive heigths with the properties h{i) = 0,1,2... and 
h{i + 1) — h{i) = +1,0,-1. This height model can also be considered as describing the 
adsorption of a polymer chain onto the wall. We further assume that the first height is fixed 
{h{i = 1) = 1). The partition function of the model reads 



/ N-l N \ 

Zn= expi-pjJ2 \h{t + l)-h{i)\+PuoY.Kfi] 

{h{i)} \ i=l 1=1 ) 



(83) 



Following closely the steps of the previous section, we first express the partition function 
Zjv as a function of its end-point {h{N) = h), and write 

Zn = Y. ZN{h) (84) 

At this stage, it is convenient to introduce some notations and relations. 

Let B{1) represent the partition function for a chain starting from h = 1 and arriving at 
/i = for the first time after I steps. We then denote by A^(/ — 1, h) the partition function 
for a chain of (/ — 1) steps going from h = 1 to h > 1 m the presence of an absorbing 
boundary at /i = 0. Setting t = e"^"^, the method of images0 yields 

B{1) = t {R{1- 1,0) -R{1- 1,2)) (85) 

and 

M{1 -l,h) = R{1 -l,h-l) - R{l~l,h + 1) (86) 

where R{1, h) denotes the partition function of the chain in the absence of the wall. We 
explicitly have 

R{1, h)=Y: Cfcf^t^ (87) 
p=\h\ 

normalized to EtL-i R{l,h) = {1 + 2t)K 

We now decompose Zj^{h) into free and adsorbed segments as 

+00 oo oo oo oo 2p 

^^W = EE E- E E^(^-(E^^ + 0) (88) 

p=0ii=l«2 = l l2p = ll=l i=l 

i3(/i)(e^""'^t)S(/3)(e^"«'^t)...(e'^"o'^^t)A^(/ - 1, h) for /i > 1 

+00 oo oo oo 2p 

p=0«i=lZ2=l hp=l i=l 
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It is again convenient to consider the grand canonical partition function 



Z{K) = ^ K^'Zn (90) 

Af=l 



Using 



+00 i_./i_4 



tK 

1-K 



2 



B{K) = = ' ' (91) 

1=1 ^l-K 

we finally have 

ZiK,2h)= , ^ ^ . \ (l^lJiBiK) - l] (92) 

1 - BiK) '^'' 7 t\ tK ^ ' ^ ' 

Z(K,2h + l)= ,^J—^—-B(K)(^—J^B(K)-l)^ (93) 

^ ' ^ 1 - B(K) ? t ^ ^ \ tK ^ ^ J ^ ^ 



B. Thermodynamics 

Equations ( P^f ) and (^ give 

Z(K) = g ft) = — ^ + (94) 

h=o 1 - B[K) ^^j^^^^^ t 2--^B{K) 

The radius of convergence -ft'*(/?) of Z{K) is either given by the solution of the equation 
(2 - ^B{K) = 0) , or by the solution of the equation (l - B{K)^^^^ = o) . The former 
corresponds to a value 

and physically describes the delocalized phase. The latter is given by 

KM = r^=^ (96) 



e'3"0-l 

The chain is localized close to the wall as long as -ft'*(/3) < Kcteioci i-e. T < = 1 / (3c with 



e--' = (97) 

1 + t ^ ^ 

m agreement with reference^!. 
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C. Density profile 

We follow equation (^6]) and define a normalized density profile by 



p{h) = lim ^^^^^ 



N 



Inverting equation (|9^, we get 



p{h) = p(0) 



2t 



in the localized phase. The characteristic length reads 



In 



1 / /e/?"0-l+4f^ 



2t \ V e/^^o-l 



and diverges at the transition as 



(98) 



(99) 



(100) 



(101) 



in agreement with equation (7.8) of referenceS. The fraction of adsorbed monomers vanishes 



as 



= 0) ^oc^/c - T 



(102) 



D. Probability distributions of adsorbed and desorbed segments 



Since the detailed procedure has been given in section [II E| , our presentation will be rather 
sketchy, and we only give the main results. Let us denote by A (resp. D) the adsorbed (resp. 
desorbed) monomers of the chain. The probability distribution of desorbed loop lengths will 
be denoted by Pd{1), and its Laplace transform by Pd{s). We get 



1=1 



(103) 



yielding, at criticality, to 



PdIs) 



T=T, ' 1 + 2t 



t 



v^ + 0(s) 



This corresponds to an algebraic decay of -Pd(0 large (/), as in equation 
in the region T — ^ Tc the moments of -Pd(0 diverge as 



(104) 
Similarly, 

(105) 
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As in the copolymer case, the correlation length along the chain scales as ( , ^ ) , in 

r- 1 \(lc-l) J 

agreement with equation (6.22) of reference". As previously found, the small loops have a 
finite weight at T^. One has for instance 

As for the probability distribution of the adsorbed segment lengths, we find an exponential 
form, even at Tc. 

Pa(/) = (1 - i^.e^"°)(ir.e^"°)'-^ (107) 



where is given in equation (p6|). 

We conclude this section on wetting by pointing out that a finite size relation for the 
specific heat, quite similar to the one derived in equation ([76|), has been found in referenced. 



IV. CONCLUSION 

We have presented a grand canonical approach to the localization of a single polymer 
chain either at a fluid-fluid interface, or at an attracting hard wall. This method, which rests 
on the absence of self avoidance in the models, provides detailed informations on the loop 
length distribution function. The extension to various disordered situationsi'S&ilwould be 
of great interest. 



16 



0.06 



0.05 



0.04 



0.03 

2, 



04 



2.06 



2.08 
T 



2.1 



2.12 



FIG. 1. Specific heat per monomer vs temperature for = 1 and qs = 0.5 in the critical 
region for 2N = 18000 (o), 14000 (□), 12000 (a), 10000 (v). The arrow denotes the critical 
temperature Tc = 2.078..., where equation (|7^ predicts C(Tc) = 0.047454.... 
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2.14 



FIG. 2. Susceptibility per monomer vs temperature for = 1 and qb = 0.5 in the critical 
region for 2N = 18000 (o), 14000 (□), 12000 (a), 10000 (v). The arrow denotes the critical 
temperature Tc = 2.078..., where equation ( p^ predicts x(^c) = 2.16289... 
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FIG. 3. Test of equation (|l]): Plot of V2N (i-^'aivCr)) vs temperature for qa = I and qb = 0.5 
for 2N = 18000 (o), 14000 (□), 12000 (a), 10000 (v). The arrow denotes the critical temperature 
Tc = 2.078..., where equation (|8l|) predicts a value 2.0279.. 
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APPENDIX A: TRANSFER MATRIX APPROACH 



In this appendix, we will briefly present the transfer matrix calculations for the periodic 
AB copolymer chain. We decompose according to the position n of the end point of 
the chain, {Zn = J2n=-oo ^N^n)), and using the Hamiltonian of equation (|1|), we obtain the 
recursion relations 

ZN+i{n) = e/^''-+i^s'^['^-^] (Z^in + 1) + Z^in - 1)) (Al) 

It is convenient to introduce the Fourier transforms for the two half-spaces 

+00 

ZUk) = T.^"'^'"'^Z^in) (A2) 



n=l 




ZN{k)= E e^'^'Z^in) (A3) 



as well as their inverse 



/-27r 

Zwin) = / — e-^^("-^)Z+(A;) for n > 1 (A4) 
JO ivr 

2^ dk 



ZN{n) = I Pe-*'='^Z;^(A;) for n<0 (A5) 
Jo 2% 

The recursion relations now read 

Zjj+^{k) = 6^"?^+^ (2cosA;Z+(A;) - e-'^Z^il) + Zjv(O)) (A6) 

Z],^i{k) = 6-""^+' {2coskZ^{k) - e^^ZMiO) - ^jv(l)) (A7) 

To study the thermodynamic limit for our problem, (g2p+i = Qa > and q2p = —qs < 0), 
we look for a stationnary solution Z^(fc) with eigenvalue A. We find 

Z+{k) = —-, ^- —[-e-"'Z{2) + e^^Z(O)] (A8) 

Z-ik) = ^^_,(,,_J_,^^^,, [(e^^^- - 1 - e-)Z(O) + e^/^-Z(2)] (A9) 



Using the inversion relations (|A5|) , we obtain the following consistency equations for Z{2) 
and Z{0) 

Z(2) = -/+(A)Z(2) + /o+(A)Z(0) 

Z{0) = e'^''-I,{X)Z{2) + {{e'^'^- - 1)1, (X) - I,{X)) Z{0) (AlO) 
whith the notations 

i}W = lJp[y = — I — j (All) 

J ( ) f^"" dk 1 1 (A12) 

° J-n 2ny-cos^k ^y{y - 1) 

^ J-n 2ny- cos2 A; \j y - I Jy(y - 1) 
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For a localized eigenstate, both Z{2) and Z{0) are non zero. As a consequence, the deter- 
minant of the system ( |A1CI|) has to vanish, and this gives an equation for the eigenvalue A 
of the localized state, which reads in our case {qa > Qb) 

where K^[j3) is given in equation (|19|). The free energy is then given by 

m = hm 4^ = -i^ (A15) 
■'^ ^ N^oo (32N 2(5 ^ ' 

in agreement with equation (^. 

As a final remark on this approach, we mention how one may find the polymer density 
profile p{2n) in the localized phase. The Fourier transforms of p{2n) on the two-half spaces 
are indeed directly related to the locahzed eigenvector given in ([A8|) by 



hoo 



V e^'^^^n-i) (2^) = . ^''^^j (A16) 

n=o Z+(0) + Z-(0) 

V e^"X2n)= . ^'^^j (A17) 

n£^oo Z+(0) + Z-(0) 

in agreement with the result (|33|) . All the results of section || can be thus recovered, except 
for the results of the loop length distribution function of section |IIE. 
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